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Abstract 

We construct the regularised Wheeler-De Witt operator demand- 
ing that the algebra of constraints of quantum gravity is anomaly free. 
We find that for a subset of all wavefunctions being integrals of scalar 
densities this condition can be satisfied. We proceed to finding ex- 
act solutions of quantum gravity being of the form of functionals of 
volume and average curvature of compact three-manifold. 
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1 Introduction 



One of the outstanding problems of modern theoretical physics is the con- 
struction of quantum theory of gravity Q. Indeed, it have been claimed 
many times that various unsolved problems like the cosmological constant 
problem, the problem of origin of the universe, the problem of black holes ra- 
diation will find their ultimate solution once this theory is finally constructed 
and properly understood. Some 0, claim that the theory of quantum gravity 
will also shed some light on the fundamental problems of quantum mechanics 
and even on the origin of mind. These all prospects are very exciting indeed, 
however, up to now, the shapes of the future theory are still very obscured. 

Nowadays there are two major ways of approaching the problem of quan- 
tum gravity. The first one is associated with the broad term 'super strings'. 
In this approach the starting point is a two-dimensional quantum field theory 
which yields quantum gravity as a part of the resulting low-energy effective 
theories. It is clear that in superstrings, like in other, less developed ap- 
proaches in whose gravity appears as an effective theory, it does not make 
sense to try to "quantize" classical gravity. 

In the canonical approach one does something opposite: the idea is to 
pick up some structures which appear already at the classical level and then 
promote them to define the quantum theory. In both the standard canoni- 
cal approach in metric representation, which we will follow here, and in the 
approach based on loop variables 0, these fundamental structures are con- 
straints of the classical canonical formalism reflecting the symmetries and 
dynamics of the theory, and their algebra. There are good reasons for such 
an approach. The equivalence principle is the main physical principle behind 
the classical theory of gravity; this principle leads to the general co-ordinate 
invariance and selects the Einstein-Hilbert action as the simplest possible 
one. 

Another building block of quantum theory is the quantization procedure. 
Here one encounters the problem as to if a generalisation of the standard 
Dirac procedure of quantization of gauge theories in hamiltonian language is 
necessary. This would be the case if one shows that the standard approach 
is not capable of producing any interesting results. It is not excluded that 
this may be eventually the result of possible failure of investigations using 
standard techniques, however, in our opinion, at the moment there is no 
reason to modify the basic principles of quantum theory. 
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Our starting point consists therefore of 

(i) The classical constraints of Einstein's gravity: the diffeomorphism con- 
straint generating diffeomorphism of the spatial three-surface "of con- 
stant time" 

V a = V b n ab (1) 
and the hamiltonian constraint generating "pushes in time direction": 

H = K 2 G abcd n ab 7v cd - ^-Vh(R + 2A) (2) 

In the formulas above 7r ab are momenta associated with the three-metric 

h a b, 

Gabcd = 2^Jl (^achbd + h ad hb c — h a bh c d) 

is the Wheeler-De Witt metric, R is the three-dimensional curvature 
scalar, k is the gravitational constant, and A the cosmological constant. 
The constraints satisfy the following Poisson bracket algebra 

[V,V]~V, (3) 

[V,-H\~H, (4) 

(5) 

(ii) The rules of quantization given by the metric representation of the 
canonical commutational relations 

7T ab (x),h cd (y)]=-i5i a 5 b j5(x,y), 

5 



n ab ( x ) = -i 



Sh ab (x) 



(iii) The Dirac procedure according to which one imposes constraints quan- 
tum mechanically by demanding that they (or, better, the correspond- 
ing operators) annihilate the subspace of the Hilbert space called the 
set of physical states. Bearing in mind the notorious regularisation and 
renormalization problems of quantum field theory, one should clearly 
state what the phrase "corresponding operators" means. In general, 
different choices of such operators could result in different quantum 
theories. 
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In the canonical approach, the points (i) to (iii) above encompass the 
whole of the input in our disposal in construction of the quantum theory. In 
particular, we do not know what is the correct physical inner product, and 
thus we do not know if the relevant operators are hermitean or not. Besides, 
we do not even know if, in the case of quantum gravity, we should demand 
these operators to be hermitean: the hamiltonian annihilates the physical 
states (the famous time problem ||) an d thus unitary evolution does not play 
any privileged role anymore. It follows that, perhaps, we cannot distinguish 
"relevant" wave functions by demanding that they are normalizable, as in 
the case of quantum mechanics, in fact, since the probabilistic interpretation 
of the "wavefunction of the universe" is doubtful, it is not clear at all if the 
norm of this wavefunction is to be 1. 

In the recent paper a class of exact solutions of the Wheeler-De Witt 
equation was found. In that paper we used the heat kernel to regularise 
the hamiltonian operator and inserted the particular operator ordering. The 
question arises what is the level of arbitrariness in this construction. In other 
words, could we construct other (possibly simpler) regularised hamiltonian 
operators and what would be their properties? This question is the subject 
of the present paper. 

It is clear from the discussion above that the only principle, we can base 
our construction on is the principle that the algebra of constraints is to be 
anomaly-free, that is, the corresponding algebra of commutators of quantum 
constraints is weakly identical with the classical one. This means that the 
structure of the Poisson bracket algebra (@-0) is to be preserved, in the sense 
which will be explained below, on the quantum level. The following section is 
devoted to the analysis of this problem. In section 3 we investigate solutions 
of the resulting equations. Some more technical results are presented in the 
Appendix. 

2 The commutator algebra and construction 
of regularised operators 

As explained in Introduction, our starting point in construction of the quan- 
tum hamiltonian operator (the Wheeler-De Witt operator) is the algebra 
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(^-§|) and we demand that the same algebra holds on the quantum level. 
At this point one should ask the question why we impose such a condition. 
One of the possible answers is that the closure of the algebra is the only 
principle which makes it possible to find operators corresponding to classical 
constraints, required by the Dirac procedure. Another argument to be found 
in the literature is that if the algebra is anomalous (that is, if there are ad- 
ditional terms resulting from the commutators of the constraint operators) 
one 

cannot find any solutions of the quantum constraints. This does not apply 
here since we know explicitly that for a particular regularisation/regularisation 
prescription introduced in || a class of solutions exist, and on solutions the 
algebra closes identically. The above argument can be therefore rephrased 
as follows. We want the algebra to close because if it does not, then the 
solutions we find will have to be in the kernel of the anomaly (we know that 
the kernel is non-empty because solutions do exist) which would mean that 
the conditions we impose will be more restrictive than the ones imported 
from the classical theory. This by itself is not a disaster, since, in any case, 
the classical limit will be the same, but we would like to depart from the 
classical theory as little as we possibly could. 

From our point of view there is another important argument in favour of 
preserving the constraint algebra structure. The vanishing of anomaly is, as 
it turns out, a quite restrictive condition which makes it possible to restrict 
the form of the employed regulator. The idea is therefore to find a regularisa- 
tion/renormalization procedure consistent with the symmetries dictated by 
the classical general relativity and to restrict it further by demanding that 
there are solutions of the theory in a class of natural wavef unctions. 

The problem of commutator algebra has been analysed in ||, |J (and 
recently in J?]]) with the result that formal manipulations involving point 
splitting lead to ambiguous final expressions. This conclusion is hardly sur- 
prising: It is well known || that to compute an anomaly one should first 
define the space of states on which the operator in question act. Then one 
should clearly state what is the procedure of extracting the finite part of 
formally divergent expressions. Thus the right question to ask is not what 
is the formal commutator of constraints but: Given a space of states, does 
it exist a regularisation/renormalization prescription such that the renor- 
malized action of the operators on the states closes? It should be stressed 
that this question is based on the basic physical interpretation of the rele- 
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vant operators; indeed, the constraints operators are generators of physical 
transformations which make sense only in terms of results of their action on 
appropriate states. As it will be seen below, the condition guaranteeing the 
absence of anomalies, not surprisingly, becomes different when the operators 
act on different states. The technical reason is simply that different states 
pick up different parts of the regulator. This fact is, of course, well known in 
investigation of anomalies in quantum field theory in canonical quantization 
language (cf. [§). 

As stressed above, we must start with choosing the initial space of wave 
functional. We assume that this space of states is the space of functions of 
Riemannian functionals, i.e., integrals over compact three-space M of scalar 
densities built of polynomials in Ricci tensor, like V = J M Vh, 1Z = J M VhR, 
etc.: 



We choose the following representation of the diffeomorphism constraint 



where we employed the notation meaning that the covariant derivative 
acts at the point x. Then we see that diffeomorphism constraint annihilates 
all the states and the commutator relation (||) is identically satisfied. This 
is the reason for a particular, natural ordering in (|7]). Moreover we see that 
the relation @ reduces to the formal relation 



Now we must turn to the heart of the problem, the construction of the 
Wheeler-De Witt operator. It is well known that second functional derivative 
acting at the same point on a local functional produces divergent result. We 
deal with this problem by making the point split in the kinetic term, to wit 



(6) 




(7) 



(8) 




dx K a b c d[X) x j i) 
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Sh ab (x) Sh c d{x') 



where K a b c d(x,x f ;t) satisfies 
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By virtue of the correspondence principle, we take 



K abcd (x, x'- 1) = G abcd (x')A(x, x'- 1) (1 + K(x, t)) , (9) 

where 

exp(-±N ab {x){x-x') a (x-x') b ) 

A & x ; *) = — 4^72 

and K(x, t) is a power series in t vanishing at t — 0. Using the fact that t 
has dimension m~ 2 we make the following expansion for K and N ab 

K{x, t) = a tR + (axR 2 + b l R oh R ab ) + . . . , (10) 

N ab (x) = h ab + 2t(A R ab + B h ab R) + t 2 (B 1 R c a R cb + A 1 R ab R + C 1 h ab R 2 ) + . . ., 

(11) 

where . . . denote the higher order terms which will not concern us, and a, b, 
A, B, C are the free parameters to be fixed. 

Next we must resolve the ordering ambiguity in the operator TC. To this 
end we add the new term L ab (x) r, 5 , r , where 

Lab = ah ab + (3h ab R + jR ab + . . . (12) 

contains free coefficients do be fixed along with the coefficients in K and N ab . 
Thus the final form of the Wheeler-De Witt operator is[] 



T-C(x) = k 2 J dx' K abcd (x, x'\ t) 



5 5 



Sh ab (x) 5h cd (x') 



+ K2i - w ^ + ^ +2A) - (13) 

To set the stage, we still need to define the action of operators on states. 
To this end we must discuss the issue of regularisation and renormalization. 
The operator (|T3|) acting on a state (defined as an integral of a scalar density) 
produces, in general, terms with arbitrary (positive and negative) powers of 
t. This provides the regularised version of the operator since all the terms 
are finite, and singularities of the form 5(0), 5'(0), etc. are traded for terms 
which are singular for t — > 0. Observe that the singular part of the action of 

1 In the paper ||] we took K a b c d(x, x') = G a bcd(x)K(x,x'), where K was a heat kernel, 
and L a b was taken to be the functional derivative of K a b c d with respect to h c d- 
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the operator on a state depends on this state. To renormalize, we follow the 
procedure proposed by Mansfield [10], based on analytic continuation, which 
result in the following: the terms with positive powers of t are dropped, and 
the singular terms of the form ^3/2 t~ k ^ 2 are replaced by the renormalization 

coefficients p k . Thus we are given a finite action of the Wheeler-De Witt 
operator on any state. 

There is a number of important comments that must be made at this 
point. It is easy to see that the singular part of the regularised action, and 
thus the renormalized action of the Wheeler-De Witt operator on a state 
does depend on the state (cf. ( |A.1| ) and ( |A.2| ) in the Appendix.) This is 
clearly a natural feature of any regularisation technique. Thus, as already 
stressed above, the condition for anomaly cancellation should be analysed 
state-by- state. 

Next, the coefficients in the regulator K are metric dependent. This 
should not be understood as an indication that the regulator depends on a 
background metric. It was observed by many that a wonderful feature of the 
metric formulation is that the metric (and its derivatives) appearing in the 
commutator is to be understood as result of action of the metric operator 
on a state. But in the metric representation the metric acts by multiplica- 
tion (h a b(x)\* >= h a b(x)\* >), and thus we can always use h a b(x) instead 
of h a b(x). Also, the regulator seemingly depends on a background structure 
through the presence of the explicit x a terms in the exponents. Such terms 
are necessarily present in any regulator based on point splitting technique. 
However it will be shown below that there is no anomaly in the quantum 
mechanical commutator of hamiltonian and diffeomorphism constraints and 
this means that the background structure dependence disappears in the final 
results. 



Now we can turn to the interpretation of equation @. According to 
our general philosophy explained above, we understand it in the following 
way. A constraint operator acts on a state and after renormalization gives 
another state depending on renormalization constants and the parameters 
of the regulator. On this resulting state the second operator acts. Thus 
the formal relation (7) is defined to mean (the state \1/ is, by definition, 
diffeomorphism-invariant) 

V{m) ren ~ (W)ren, (14) 
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and, similarly, for the hamiltonian-hamiltonian commutator 

(H[N] {H[M}m) ren ) ren - (Ti[M] (H[N}*) ren ) ren = (15) 

for all M and N. In the formula above we used the smeared form of the 
Wheeler-De Witt operator 

H[M] = J dxM(x)H(x). 

Let us turn back to equation (|14|). Since the action of diffeomorphism 
is standard, it suffices to check that {7i^) ren is a scalar density. But this is 
clearly the case: the first functional derivative acting on a state produces a 
tensor density T ab (x'). After acting by the second derivative and contracting 
indices, we obtain the terms of the form 

T (x')8(x', x)+T 1 (x')oX7 x 'oV x '5(x', x)+T 2 {x')oV x ' oV 1 ' oV 1 ' oV*' ' 5{x' , x)+... 

where o denotes various indices contractions, and T n are tensor densities. 
These terms are multiplied by A(x, x'\ t) and integrated over x' . Now we in- 
tegrate by parts which results in replacing covariant derivatives acting on K 
with appropriate powers of t~ x multiplied by some coefficients. After renor- 
malization we obtain a scalar density as required. The action of the L term 
clearly gives the same result. Thus 

For the states being integrals of scalar densities there is no anomaly in 
the diffeomorphism — hamiltonian commutator 

This result is quite important because the anomaly in the string theory 
appears in the diffeomorphism — hamiltonian commutator. It proves also 
that, in spite of implicit co-ordinate system present in the construction of the 
regulator, the three dimensional diffeomorphisms are not broken by quantum 
corrections. 



Now we turn to the most complicated problem, the hamiltonian — hamil- 
tonian commutator (15). Our goal will be to use this equation to partially 
fix the free coefficients in K, N a i, L a ^. These coefficients will be further fixed 
by demanding existence of solutions of Wheeler-De Witt equation. In what 
follows we will be interested in solutions of the form ty(V,TZ). Therefore we 
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check explicitly the closure of the algebra only for the states of this form. 
We will comment on the general case at the end of this section. 

Let us start with the simplest state \? = 1. The action of the first smeared 
operator gives simply 

(H\M]y) ren = — / dx Vh(x)M(x)(R(x) + 2A). (16) 

K J 

Now we have to act on the right hand side of the above equation with the 
operator (TC[N], then renormalize the result, and finally subtract the result 
of the same calculation with N interchanged with M. After rather tedious 
computation one finds in the commutator the term proportional to [N, M] a = 
NV a M - MV a M which must vanish, to wit 

p w a^-V a R - (V a L - V b L b a ) = 0, (17) 

where L = h ab L ab . Using Bianchi identity and the expansion ([TJD we find the 
first relation between coefficients, to wit (it will soon turn out that . . . terms 
in L ab vanish) 

^(3p«ao- 7 )-2/3 = 0. (18) 

Now let us turn to the states depending of V = f M d 3 xyh. Let H[M] 
act on this state. From Eq. ( |A.1| ) we see that we have an equation which is 
of the form flT6|) (with different coefficients which will include V.) Therefore 
the condition for vanishing of the commutator is the same as above, (|Hf). 

Let us pause for a moment with investigation of the algebra to make an 
important observation. We want some \I/(V) to be a solution of the Wheeler- 
De Witt equation. From the computations above we see that the double 
derivative term in 7i will produce terms up to order R. It follows that, while 
solving the equation, we would not be able to cancel terms of higher order 
in R (like R 2 .) Therefore, all the terms in L ab expansion (|T2|) denoted by . . . 
must vanish. Thus we take 

L ab = ah ab + (^R ab + (3h ab R), 

where the coefficients (3 and 7 are subject to the condition (|18D. 

Now we turn to the wavefunction \& = $?(1Z). Let us analyse the action of 
the commutator of hamiltonian in a number of steps. The first observation 
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follows from the term in (|A.2 ). It can be checked that after acting on 



this term by 7i[N] one obtains a term proportional to W which contains 
unremovable anomaly of the form 

3 3 



[N, M} a (-R ab V b R - Yq RV ^ R ) 



Till this point we assumed that the wavefunction ^/(TZ) was arbitrary, thus 
anomaly multiplying W was to vanish independently of possible anomalies 
multiplying different derivatives of But this is, clearly, cannot be accom- 
plished. We thus have no choice but to restrict It would seem that fixing 
the background geometry 

- \r 2 + R ab R ab = (19) 

o 

would do, but this cannot be done because some external condition may be 
applied only after the commutator is computed, and not at the first step. 
The only way out is to make ^" = 0, or proportional to (where the 
terms quadratic in curvature are already present.) This means that either 
$?(1Z) = ATI or ty(TZ) is a linear combinations of exponents exp(ulZ). It 
should be stressed that since solutions we are after will necessarily have the 
form of exponents, the restriction we are making is not as severe as it would 
seem at the first sightQ. The inspection of equation ( |A.2| ) clearly shows that 
the most economic way is to take B = —\J- In this way we can cancel the 
anomaly for arbitrary u. 

It turns out that the anomaly is proportional to [N, M] a times a combi- 
nation of five different tensorial objects. To cancel the anomaly proportional 
to R b c VbRa we must put B\ = 0. Similarly, the condition for the anomaly 
proportional to i^V^i?^ to vanish is b\ = 0. From the conditions for R^VbR 
and (Va-R)-R, and the equation B = — |j7" we find expressions for Ai, C\, and 
ai, to wit 

11 1 /5 



A 1 = -a A + — a + (^-7 + (3) , (20) 

9 A 29 7 
C\ = -~a A - — a - - (21) 



16 
1 

66~ u n~ u " u ' 33 



29 9 , 1 f\A n 37 . 
o-i = ~7^ a o ~ TT a oA) ttt 77P + 77777 (22) 



2 It is sufficient to check that there is no anomaly in an neighbourhood of solutions. 
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There is one equation remaining, being a coefficient of V a -R anomaly which 
relates o parameters to each other: 

^(-s-i*-T* + 5*H a = a (23) 

We are left therefore with six free coefficients of the regularised Wheeler-De 
Witt operator A , B , a , a, j3, and 7 subject to two linear equations (|I8|) 
and (^). 

Thus the final form of the regularised Wheeler-De Witt operator which 
preserves the constraint algebra is (to linear order in R and with included) 

H(x) = k 2 [ dx' G abcd (x') A(x,x';t)(l + a tR + ...) — - -— — + 

5 1 

+ k 2 (ah ab + {(3h ab R + jR^)) (x) + — Vh(R + 2A). (24) 

The formula ([24]) completes our construction of the Wheeler-De Witt 
operator. As compared to the choice made in the paper 0, where we used 
the heat kernel and L ab was its functional derivative, here we gained much 
more freedom in the form of additional free constants. These constant will be 
further fixed by demanding that the Wheeler-De Witt equations possesses a 
maximal number of solutions, that is that there are solutions ^(V), ^(TZ), 
and #(V,ft). 

It is possible to extend the above analysis to the states being functionals 
of higher powers of curvatures. To this end one has to add terms of order 
t 3 to the regulator K, compute the commutator, and fix the coefficients as 
it was done above. It seems quite likely that the resulting equations could 
be solved. However, the computations are becoming extremely tedious, and 
for that reason in this paper we were not able to address the question of 
anomalies for higher states. 



3 Solutions 

From the previous section we know that the most general form of the Wheeler- 
De Witt operator satisfying our criteria is given by equation (0). Now, em- 
ploying this operator, we will try to find a class of solutions of the Wheeler-De 
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Witt equation. It should be stressed at this point that we regard the exis- 
tence of a maximal possible space of solutions as an ultimate condition fixing 
the operator completely. The reason for is that any regulator defines a quan- 
tum theory with a certain set of solutions. Clearly, a theory with richest set 
of solutions is most interesting. Thus our goal is twofold: to find solutions 
and to fix the operator as to allow for the maximal possible number of them. 

We will consider the states of the form \P = ^(V, TV). It is clear from the 
form of the Wheeler-De Witt equation that the resulting equations, multiply- 
ing various scalar functions will be linear, therefore, without loss of generality, 
we can assume that ^(V, 1Z) = exp(aV + ujIZ). 

First we solve the equation for the coefficient u. From the part of the 
equation involving the square curvature terms, ( |A.2| ), we see easily that since 



B = —\J, w = — J = 7. Let us then turn to the coefficients multiplying 



k 2 (aX + uy- -toa^j + = 0, (25) 
y = p^ ( Ao + 3B -^a + 7 X ' 



where 



, OL. 

2 87 2 
The coefficient multiplying the y/h term reads 

K * f_3 a 2 _ 21^(3) + l aa _ 3 p(5) \ + 2A = Q _ (2g) 

2 2 / k 



3 Observe that in addition to terms presented in appendix, ( |A.2[ ), ( |A.l| ), there is another 
term resulting from the action of second functional derivative, one on V, and one on 1Z. 
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Solutions with A ^ 
We consider three cases: 



Case I. u> = 0. In this case we have 

^ = exp(5"V), a 



K 4 X ' 



and equation (26) gives a condition for the parameters of the regulator which 
can be solved for a. 

Case II. oj — 7 and a = 0. We have 

1 

y- 



k 4 7 



and the condition relating 7 to the bare coupling and renormalization con- 
stants 

1 = 5?JW- (27) 

Case I/Iu; = 7 and a 7^ 0. We find a = 4a — 7// 3 ) and the condition X = 
27. This condition can be solved along with conditions from the previous 
section to give expressions for the regulator parameters. 

Collecting all results we finally have 

j = exp f- 3 ^-V^ Case I; (28) 

9 n = exp (g^j-fcj Case II; (29) 

/ 4A \ 

^ = ex p(- 3 X V+ 3^)^J CaSeIIL (30) 

Observe that solution III is a product of solutions I and II. We will return to 
this observation below. 

Solutions with A = 
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it is easy to see that (assuming p^ ^ 0) in this case the wavefunction 
^ = exp(jlZ) does not solve the Wheeler-De Witt equation. Taking X = ^7 
as above we find the solution 

= exp (-ijv) (31) 

The second solution is of the form 

V IIIo = exp (<xV + l1Z) , (32) 
where a is a solution of the following quadratic equation 

a 2 + -^-a - 4p (5) 7 = 0. (33) 

Depending on the value of 




we have either two real, or two complex, or one real solution. 

Thus we have three different wavefunctions (for both cases A = and 
A ^ 0) being solutions of the Wheeler-De Witt equation and containing 
functionals of order at most linear in R. It is very interesting that the 
solutions depend on the bare coupling constants k and A and only on a 
single renormalization constant p^ . Of course, any linear combination (with 
complex coefficients) of the solution is a solution. Such combinations will be 
called below "Schrodinger cat universes". It can be argued that, contrary to 
the real solutions, complex solutions will in general possess a nontrivial time 
evolution. 



A Renormalized action of H[M] 

Here we present the calculation of the renormalized action of hamiltonian 
constraint on states. We have 

5 S 21 

GabcdW) v = -—K x - O- 

5h ab (x) 5h cd {x') 8 
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Thus 

/ rfV K(x - x '- W^ hhlwsh!^ = 
-^T WW {w + W'"" R + 0(f) ) • 

Using this result and renormalizing, we obtain 

H[M]V(V) = J d 3 xVhM^-K 2 -^"(V) + 

+^'{V) (-y (P (3) + aoP^R) + \l)+±- 2 {R + 2A)*(V)} , (A.l) 

where L = L ab h ab . 
Similarly, 

G abcd (x , )—^ T1 — ^— -TZ = ? -R(x')S(x - x') + U x ,5{x - x 1 ). 
bh ab {x) bh cd {x') 8 



Thus 



/ rfV K(x,x';t)G abcd (x') K 
J ohnhix) bhrdix') 



Vh 



5h ab (x) 5h cd (x') 

3 



2 

'7 



+ p(3)i^o + 35o-^ao + £) + 



+p« |i? 2 (-ao + A 1 + 3d - -ai + (A + 350)^) + 
+ R ab R ab (5!-^6i + 3£)i)} 



Thus we obtain 



+K 2 tt'(ft) \BR 2 + JR ab R ab + R (p (3) (Ao + 35 - \a Q + £) + ±a) - 

±(R + 2AWJl)\, (A.2) 



+ 
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where 



B 



-ai + (A + 3B )a + A 1 + 3C ± 



and 



J 



7- 



(A.3) 



(A.4) 



Equations (|A.1|) and ( |A.2j ) are basic for our investigations in the main body 
of the paper. The expressions in the parentheses { * } in these equations are 
the Wheeler-De Witt equations for the corresponding wavef unctions. 
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